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I. Inleiding 


Statistische fysica houdt zich bezig met de fysische eigenschappen van 
macroscopische systemen, zoals een stuk metaal, een cylinder met gas, of 
een vat met vloeistof. Het essentiële kenmerk van macroscopische systemen 
is dat het aantal deeltjes (atomen, moleculen) zeer groot is. Statistische 
methoden worden beter bruikbaar naarmate er meer deeltjes zijn, en naarmate 


we minder geinteresseerd zijn in het specifieke gedrag van een enkel indi- 


vidueel molecuul. In de praktijk is het deeltjesaantal meestal van de orde 


van 10° of meer. De toestand en de eigenschappen van een macroscopisch 

systeem worden gekarakteriseerd door een gering aantal macroscopische grootheden, 
zoals de temperatuur, de druk, de electrische stroom, het geleidingsvermogen, 

de thermische uitzettingscoefficient, de warmtecapaciteit. De statistische 

fysica probeert nu deze grootheden af be Teile uit de eigenschappen van een of 


enkele moleculen, bij voorbeeld de massa of het electrisch dipoolmoment van 


een molecuul, of de wisselwerkingspotentiaal tussen twee of meer moleculen. 


Niet altijd is het duidelijk hoe een macroscopische grootheid in termen van 
moleculaire eigenschappen gedefinieerd zou moeten worden. 

Het is niet moeilijk te begrijpen waarom statistische methoden nodig zijn. 
Als we willen onderzoeken hoe een macroscopische grootheid verandert als functie 
van de tijd, zouden we kunnen proberen om mechanica toe te passen. Als we de 
wisselwerkingspotentiaal tussen de moleculen kennen, kunnen we de bewegings= 
vergelijkingen voor alle deeltjes opschrijven (de wet van Newton als we klassieke 
mechanica gebruiken, of de Schrödinger=vergelijking als we de quantum=mechanica 
willen toepassen). Deze bewegingsvergelijkingen moeten we dan oplossen bij een 
gegeven beginsituatie. Het is duidelijk dat deze methode niet uitvoerbaar Is 
als het deeltjesaantal van de orde van 10°? is. Alle computers ter wereld zouden 
met elkaar de beginvoorwaarden niet kunnen opnemen in hun geheugen. Daar komt 
nog bij dat we helemaal niet geïnteresseerd zijn In de precieze beweging van 
alle deeltjes. Want een macroscopische grootheid wordt bepaald door het gemiddelde 
gedrag van vele deeltjes. Als we de druk op een wand van een cylinder meten, dan 
moeten we de impulsoverdracht gedurende een korte tijd, waarin vele deeltjes tegen 
de wand botsen. In feite meten we slechts een gemiddeld gedrag. Met heel gevoelige 
apparatuur is het weliswaar mogelijk om fluctuaties in het gemiddelde gedrag waar 
te nemen. Maar ook die fluctuaties kunnen we beschrijven met behulp van gemiddel- 
den over deeltjes. Zo kunnen we begrijpen waarom statistische methoden zo belang- 
rijk zijn in de fysica van macroscopische systemen. | 

De oudste toepassing van statistische methoden in de fysica van macroscopische 
systemen was in het gebied van de thermodynamica, zoals de relatie tussen de 


temperatuur en de druk van een gas. 
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II. De eerste hoofdwet en de energie. zi 

We beschouwen een macroscopisch systeem, b.v. een gas in een cylinder, 
of een vat met vloeistof. De totale energie E van dit systeem is constant 
zolang er van buitenaf geen energie wordt toegevoerd of afgevoerd. Deze 
totale energie is de som van de kinetische en potentiële energie van alle 
afzonderlijke deeltjes. Er zijn verschillende manieren om energie aan het 
systeem toe te voeren. De meest voor de hand. liggende manier is om arbeid 
op het systeem uit te oefenen, bv. door erin te roeren, of door een kracht 
uit te oefenen op een zuiger die het systeem van de buitenwereld afsluit, 
of door middel van uitwendige electrische of magnetische velden. Maar ar= 
beid is niet de enige manier om de energie van het systeem te veranderen. 
Als de toestand van het systeem verandert, dan is de energlieverandering 
AE niet altijd gelijk aan de op het systeem verrichte arbeid W. We de- 


finieren nu de resterende energieverandering als de toegevoerde warmte Q 
Q= AE = W ae) 


Vanzelf heeft de grootheid toegevoerde warmte dan de dimensie van energie. 
De meest gebruikelijke manier om warmte toe te voeren (d.i. om energie 
toe te voeren anders dan door arbeid) is om het systeem in contact te 
brengen met een ander systeem van hogere temperatuur, bv. een gasvlam. 
Overigens hebben we het begrip temperatuur nog niet gedefinieerd. We ge- 
bruiken hier slechts het ervaringsfeit dat energietoevoer mogelijk is op 
andere wijze dan door mechanische arbeid. 

Een ander ervaringsfeit is dat het mogelijk is om een systeem thermisch 
te isoleren, door het te omgeven met wanden die geen warmte doorlaten 
(adiabatische wanden). Voor zo'n geisoleerd systeem geldt Q = O0, en der= 
halve AE = W. Dat is het geval in de klassieke calorimeter, zolang die 
tenminste gesloten blijft. De energieveranderingen W en Q in (2.1) zijn ma- 
croscopische begrippen: W is de energieverandering door op het systeem 
verrichte arbeid met macroscopische middelen, zoals zuigers, hamers, roer= 
staven, of macroscopische uitwendige velden. Op moleculair niveau is geen 
wezenlijk onderscheid tussen verrichte arbeid en toegevoerde warmte: 
beide vormen van energieverandering wijzigen de snelheden of posities van 
de moleculen. Vergelijking (2.1) definieert in onze redenering het begrip 
toegevoerde warmte. Historisch staat (2.1) bekend als de Eerste Hoofdwet 
van de thermodynamica, die behoud van energie uitdrukt in termen van 


twee meetbare, macroscopisch verschillende vormen van energieverandering 
Qen W. 
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De totale energie E ven een systeem is een grootheid die bepaald is door 
de toestand van het systeem onafhankelijk van de wijze wenken ket systeem 
in die toestand gekomen is. Een dergelijke grootheid noemen we een toe= 
standsgrootheid. De toegevoerde warmte Q en de op het systeem verrichte 
arbeid W zijn geen toestandsgrootheden. We kunnen aan de toestand van een 
systeem niet zien of die toestand (uitgaande van een zekere begintoestand) 
bereikt is door warmte toe te voeren of door arbeid op het systeem te 
verrichten. Daarom kunnen we ook niet spreken van de totale hoeveelheid 
warmte of de totale hoeveelheid arbeid van een toestand van het systeem. 

De Eerste Hoofdwet voor een infinitesimeal kleine verandering wordt ge= 


schreven als 
dE = &Q + &W | ; (2,2) 


Het symbool & geeft aan dat de betreffende verandering geen verandering in 
een toestandsgrootheid is. 

Een zeer gebruikelijke manier om arbeid op een gasvormig systeem te 
verrichten is door het volume te veranderen. Dat kan door het systeem 
in een cylinder af te sluiten met een beweegbare zuiger, zoals voorkomt 
in verbrandingsmotoren en stoommachines. Het systeem kan omgeven zijn 
door adiabatische wanden. Bij elke verandering geldt dan @Q= o, zodat 
volgens (2,2) dE = &W. We spreken dan van een adiabatische verandering. 
Anders is het als we het systeem in warmte-contact brengen met een veel 
groter systeem bij een bepaalde temperatuur. Na elke verandering zal dan 
zodanige warmteuitwisseling tussen het systeem en het grotere systeem 
plaatsvinden dat de temperatuur ongewijzigd blijft. Een verandering 
waarbij de temperatuur niet verandert noemen we isotherm. Het grote 


systeem heet warmtebad. 
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Als ik op de zuiger druk met een druk p = 0 (K is de kracht, O het opper= 
vlak), en als onder invloed van die kracht het volume V verandert met een 


waarde av, aan is de op het systeem verrichte arbeid 
aw = -p dv, | (2.3) 


Bij compressie is dV negatief, bij expansie positief. Toestandsveranderingen 
kunnen we bewerkstelligen met behulp van gewichten ‘op de zuiger, of via 
warmtecontact met een systeem van hogere of lagere temperaturen. Als we 
een systeem na een toestandsverandering aan zichzelf overlaten, dan zal 
na enige tijd de toestand niet meer macroscopisch waarneembaar veranderen. 
We spreken dan van thermisch evenwicht. We zullen nog zien dat in thermisch 
evenwicht (of kortweg: in evenwicht) de temperatuur overal in het systeem 
gelijk is. Een toestandsverandering noemen we quasi=statisch als die z0 
langzaam verloopt dat op elk moment het systeem in evenwicht is. Een quasi= 
statische verandering is dus een continue opeenvolging van evenwichtstoe= 
standen. Als we het gas comprimeren door plotseling een zwaar gewicht op 
de zuiger te laten vallen en vervolgens het systeem tot evenwicht laten 
komen, dan is er tijdens het proces geen evenwicht, en de toestandsveran= 
dering is niet quasistatisch. Alleen als we de kracht op de zuiger lang= 
zaam en in kleine stapjes laten toenemen, dan zal op elk moment het systeem 
vrijwel in evenwicht zijn. Als we nu de kracht steeds langzamer laten | 
toenemen, dan wordt een quasistatisch proces steeds beter benaderd. Wanneer 
de zuiger wrijvingsloos beweegt, dan is de druk die we op een zeker moment 
tijdens quasistatische compressie moeten toepassen gelijk aan de gasdruk op 
dat moment. Immers, als we toestaan dat de zuiger willekeurig langzaam 
beweegt, dan is een willekeurige kleine overdruk voldoende. 

Een toestandsverandering die quasistatisch verloopt en zonder wrijvings= 
effecten: noemen we reversibel of omkeerbaar. Bij een infinitesimale 


reversibele volumeverandering geldt dus 
ÈW = =pdV (ak) 


waar p de druk van het gas is. Een reversibel proces kan door een in= 
finitesimale verandering in de uitwendige krachten (in ons geval de kracht 
op de zuiger) in omgekeerde richting verlopen. Immers, de infinitesimaal 
kleine overdruk behoeft slechts veranderd te worden in een infinitesimaal 
kleine onderdruk om quasistatische reversibele compressie te veranderen | 
in quasistatische reversibele expandie. Een proces dat niet quasistatisch 
verloopt, of dat met wrijving gepaard gaat, kan niet in omgekeerde richting 


verlopen door een infinitesimale verandering in de uitwendige krachten. 
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Bij zo'n proces moet de gelijkheid (2.b) vervangen worden door de ongelijkheid 
dW > —pdV (2.5) 


Want bij niet-quasistatische compressie moet ik een druk P, Op de zuiger 
toepassen die groter is dan de gasdruk p‚, teneinde het proces met eindige 
snelheid te laten verlopen. Bij wrijving moet de druk op de zuiger Pp, groter 
Zijn dan P, teneinde de wrijvingskracht te overwinnen. Omdat P> P geldt 
dW = -P‚dV > -pdV (dV is negatief). 
Bij niet-quasistatische expansie moet Po echt kleiner zijn dan p . Voor: de (nu 
negatieve) arbeid @W geldt dan weer dW= -P‚dV > -pdV. Bij irreversibele compressie 
moet meerarbeid aan het systeem toegevoerd worden dan er weer vrijkomt bij daar- 
op volgende irreversibele expansie. 
Ruw gezegd, een irreversibel proces leidt tot verlies, van arbeid. Daaruit volgt 
direct de bovengenoemde bewering dat een irreversibel proces niet in omgekeerde 
richting kan verlopen door alleen een infinitesimale verandering in de krachten. 
De ongelijkheid (2.5) geldt algemeen voor irreversibele processen. Reële pro-= 
cessen zijn altijd irreversibel. Reversibele processen kunnen willekeurig dieht 
benaderd worden door heel kleine uitwendige overdrukken toe te passen. Voor rever- 


sibele processen kunnen we de Eerste Hoofdwet schrijven als 


dE = 8Q - pdV (2.6) 


Voor irreversibele processen geldt valgens (2.5) 


dE > dq - pav (2.7) 
De extra arbeid door overdruk die een uitwendige kracht verricht bij compressie 
draagt dus bij tot dE. 
Bij de beschouwingen in dit hoofdstuk namen we aan dat arbeid op het systeem 
verricht kan worden door volumeverandering. Als andere arbeidverrichtende | 
mechanismen (zoals uitwendige velden) een rol spelen, leiden soortgelijke beschou- 


wingen tot dezelfde conclusies. 


III. De tweede hoofdwet en de entropie 


In het vorige hoofdstuk noemden we het begrip thermisch evenwicht. In een sys- 
teem waarin de temperatuur van plaats tot plaats verschilt, of waarin de deel- 
tjesdichtheid of de druk grillig variëert als functie van de plaats heerst geen 
evenwicht. Er zal een warmtestroom ontstaan van hoge temperatuur; en een deeltjes- 
stroom, in het algemeen van hogere dichtheid naar lagere dichtheid (eventueel 
beïnvloed door uitwendige krachten), totdat na korte tijd de situatie niet meer 


verandert. Dan spreken we van thermisch evenwicht. Dit is een tYP isen macrosco 





or 


pisch begrip. Op moleculaire schaal gebeurt er nog wel degelijk lets, en de 
moleculen bewegen net zo grillig door elkaar als tevoren. De moleculaire 
(of mieroscopische) toestand van het systeem verandert nog voortdurend, en die 
wordt ook niet stationair. Het is duidelijk dat we met de macroscopische toe- 
stand (kortweg macrotoestand) en de microscopische toestand (kortweg miecrotoestand) 
iets geheel verschillends bedoelen. Een macrotoestand is gespecificeerd door een 
beperkt aantal macroscopische grootheden, zoals drukverdeling, dichtheidsverde- 
ling, temperatuurverdeling, deeltjesstroom. Het hangt af van de nauwkeurigheid 
van de instrumenten, een hoe nauwkeurige specificatie van een macroscopische 
toestand zinvol is. De definitie van het begrip macrotoestand is kennelijk wat 
vaag. Dat geldt niet voor het begrip microtoestand, gedefiniëerd als een volledige 
beschrijving van de moleculaire toestand van het systeem. In de klassieke mechanica 
betekent dit een precieze opgave van de posities en snelheden van alle deeltjes, 
bij gebruik van quantummechanica wordt de miecrotoestand gegeven door de golffunctie 
van het gehele systeem, Bij een gegeven macrotoestand behoren nog geweldig vele 
microtoestanden. Dit idee van een macrotoestand waarmee vele microtoestanden 
corresponderen is gebruikt door Boltzmann om het begrip entropie te definiëren 
in termen van de moleculaire eigenschappen van het systeem. De entropie was eerder 
ingevoerd door Clausius met behulp van alleen macroscopische grootheden. 
Boltzmanns definitie is veel meer aanschouwelijk, en we zullen ons vooral tot 
deze beperken. 
Beschouw een geheel afgesloten systeem, dus met vaste, adiabatische wanden. 
Dan is er geen energieuitwisseling met de omgeving mogelijk. Het deeltjesaantal 
N, het volume V en de totale energie E zijn dan constant. Als het systeem in 
thermisch evenwicht is, bepalen deze drie grootheden E‚N,V de macrotoestand vol- 
komen. Maar er zijn ook macrotoestanden waarbij het systeem niet in evenwicht iS. 
Dan moeten nog meer parameters gespecificeerd worden, die we kortweg aanduiden 
met ad. Een macrotoestand wordt dus aangeduid met (E‚N,V‚a). Met een willekeurige 
macrotoestand (E‚N,V‚a) komen nu vele microtoestanden overeen, waar een micro- 
toestand een eigentoestand is van de totale moleculaire Hamiltoniean. We voeren 
nu in het statistisch gewicht Q (E‚N,V‚a) als het aantal microtoestanden dat 
correspondeert met de macrotoestand (E‚N,V,a). 
Boltzmann voerde twee postulaten in betreffende het statistisch gewicht van 
een macrotoestand, die het mogelijk maken waarschijnlijkheidsverdelingen over 
microtoestanden af te leiden. We zullen beide postulaten neen en enigszins 


toelichten. 


La 


Het eerste postulaat van Boltzmann luidt: 


Voor een geïsoleerd systeem bij vaste gegeven waarden van E‚N‚,V zijn alle 
microtoestanden die deze waarden hebben even waarschijnlijk. 


pige postulaat houdt een aanname in voor de waarschijnlijkheidsverdeling over de 
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mogelijke mierotoestanden. Doordat de energie een gespecificeerde waarde heeft , 

Is slechts een beperkte verzameling van microtoestanden mogelijk, en wel de 
microtoestanden met deze bepaalde energie. De overige microtoestanden zijn 
uitgesloten. Boltzmann neemt dus aan dat de waarschijnlijkheid eerlijk verdeeld 

is over de mogelijke microtoestanden. 

Een onmiddelijk gevolg van dit postulaat is dat de waarschijnlijkheid van een 
zekere macrotoestand (E‚N,V,a) bij de gegeven waarden van E‚N,V, en voor varia- 
bele a evenredig is met het statistisch gewicht Q(E,N,V‚a) van die macrotoestand. 
In dat eerste postulaat van Boltzmann wordt het begrip waarschijnlijkheid voor 

het eerst genoemd, Het idee is dat onze macroscopische kennis over de tostas van 
het systeem niet voldoende is om ook de precieze microscopische toestand te kennen. 
Over die microtoestand kunnen we slechts spreken in termen van waarschijnlijkheid. 
Bovendien kunnen in een geïsoleerd systeem in de loop van de tijd vele verschil- 
lende macrotoestanden voorkomen, getypeerd door de mogelijke waarden van a bij 

de vaste waarde van E‚N,V. Boltzmanns eerste postulaat zegt ook iets over de 
waarschijnlijkheid om bij steekproefsgewijze metingen in zo'n geïsoleerd systeem 
een bepaalde waarde van a te meten, dus om een bepaalde macrotoestand te vinden. 
Het tweede postulaat van Boltzmann luidt: 


De evenwichtstoestand van een geïsoleerd systeem is de macrotoestand waar- 


voor het statistisch gewicht Q(E,V‚N,a) bij vaste E‚V.N. maximaal is. 


Volgens het eerste postulaat is dan de evenwichtstoestand van een geïsoleerd 
systeem de meest waarschijnlijke macrotoestand. Dat ligt wel enigszins voor de 
hand. Bij macroscopische systemen met hun geweldig grote aantallen deeltjes 
is het maximale statistisch gewicht ook geweldig veel groter dan de andere statis- 
tische gewichten. De kans op een relatief zeer kleine fluctuatie uit de macrotoes- 
tand met maximale Q is dan al miniem, en dat is juist wat we van de evenwichts- 
toestand verwachten. De evenwichtstoestand is dus niet volkomen stationair. 
Heel kleine fluctuaties komen voor, maar die fluctuaties worden relatief kleiner 
naarmate N groter is. 

We willen deze postulaten, en de begrippen microtoestand en macrotoestand 
illustreren aan de hand van het voorbeeld van een klassiek ideaal gas, dat is 
een gas van niet-wisselwerkende deeltjes. De deeltjes hebben massa m. Ze kunnen 
botsen met de wand, en ze bewegen onder invloed van een uitwendig veld (bv een 
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waar vr, ) de potentiele energie Is die het uitwendige veld beschrijft. De impuls 
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en de positie van een deeltje kunnen we aangeven als een vector (p‚r) in een 6- 
dimensionale ruimte (de zogenaamde u-ruimte). We verdelen nu deze 6-dimensionale 


ruimte in gelijke hokjes, die zo klein zijn dat de plaats en impuls van een deeltje 
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voldoende nauwkeurig gespecificeerd is door te zeggen in welk hokje de 

6-vector (D‚r) zit, maar toeh zo groot dat in het relevante gebied van de 
u-ruimte vele deeltjes in hetzelfde hokje horen. De hokjes nummeren we met 

de index a = 1,2,..:. De mierotoestand van het systeem kennen we als we weten 

in welk hokje deeltje 1 zit, In welk hokje deeltje 2 zit, enz. Laten we nu 
afspreken dat we de macrotoestand karakteriseren door te zeggen hoeveel deeltjes 
er in hokje 1 zitten, hoeveel er in hokje 2 zitten, enz. 

Dus de macrotoestand speciferen we door aan elk hokje dad een aantal deeltjes 

n toe te kennen. Dus de getallen Ín_} met elkaar spelen de rol van de variabelen 
a. Ze moeten voldoen aan twee voorwaarden, namelijk dat het totale deeltjes- 


aantal N, en de totale energie E Is. In formule 


En =N | | 13.2) 
a a 

Ene =E, | (3.3) 
a a a 


Het statistische gewicht van de macrotoestand Ed is nu het aantal manieren 


waarop ik de N deeltjes kan verdelen in een groep van N, deeltjes, een groep 


1 
van N, deeltjes, enz. Dit is de multinomiaalcoefficient 


(in, }) Epe (3) 
Volgens het tweede postulaat is de evenwichtstoestand de macrotoestand met 
maximale Q. Dus we moeten (3.h) maximaal maken, onder de bijvoorwaarden 
(3,2) en (3.3). Evengoed kunnen we log Q maximaliseren. We doen dit met 
behulp van de multiplicatorenmethode van Lagrange. Omdat N, Nys Doses: 
grote getallen zijn, mogen we log n! benaderen met de Stirling-formule 


log n! © n log n -n. (3.5) 
| continue | 
Bovendien kunnen we de n, wel behandelen alsof het variabelen zijn. Volgens 


de Lagrange-methode vinden we het maximum van log Q met de bijvoorwaarden 


(3.2) en (3.3) door te schrijven 


al log 2 - ME n, - N-B(Ze n, -E) J =o (36) 
a a 


Hier zijn \ en Bg de Lagrange-multiplicatoren. De waarde van À en B wordt 
bepaald door de bijvoorwaarden (3,2) en (3.3). Door (3.h) en (3.5) in te 


vullen in (3.6) vinden we e 
- z(log Dkt Be )dn,= o | te Sat) 


voor elke variatie dn,» Dus 
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De waarden van A en B kunnen we vinden uit (3.2) en (3.3). Voor ons is het 
voldoende te constateren dat in evenwicht de dichtheid in de u-ruimte 
evenredig is met exp(-Be). Dus het aantal deeltjes met impuls tussen 


> > > en 4 > > > 
pen p + dp en met positie tussen r en r + dr 1s 


>> 


f(D‚r)apdr = C exp -B(B /2m + v(r))] apar (3.9) 

De normering C is bepaald door de eis 
EE de ie 

(e(B,Paber = N. (3.10) 


De functie f is de Maxwell-Boltzmannverdeling. Dit voorbeeld laat zien dat 
we met de postulaten van Boltzmann de waarschijnlijkheidsverdeling voor de 
evenwichtstoestand van een geisoleerd systeem kunnen vinden. Het invoeren 
van de hokjes in de p-ruimte is een truc die nodig is om het aantal 
microtoestanden aftelbaar te maken. Bij gebruik van quantummechanica is 
dat niet nodig. 

Het is in het algemeen praktischer om met log Q te werken dan met @ 


zelf. Daarom voeren we In de grootheid S, door de gelijkheid 
S = k log @ ED 


waar k een constante is, die de constante van Boltzmann genoemd wordt. 
Later zullen we het verband leggen tussen k en meetbare grootheden. De 
grootheid S noemen we de entropie. De entropie is dus groter naarmate het 
statistisch gewicht groter is. Het tweede postulaat kunnen we nu ook 
formuleren door te zeggen: de evenwichtstoestand van een geisoleerd systeem 
Is de toestand met maximale entropie. Merk op dat S net als Q een functie 
is van (E,V‚N,a). 

Als op een zeker begintijdstip het systeem niet in evenwicht is, ver- 
keert het in een macrotoestand met niet-maximaal statistisch gewicht. 

Dan is het het meest waarschijnlijk dat het spoedig zo zal veranderen dat 
het terecht komt in een meer waarschijnlijke toestand, d.i. een toestand met 
grotere Q. Het is aannemelijk dat het statistisch gewicht van een systeem 
dat niet in evenwicht is, alleen maar groter zal worden. Dit wordt uit- 


gedrukt door de Tweede Hoofdwet van de thermodynamica: 


De entropie van een systeem dat niet in evenwicht 1s, neemt altijd 
toe. In de evenwichtstoestand bereikt de entropie zijn maximale waarde. 


De Tweede Hoofdwet zegt ongeveer dat bij een reëel proces de toestand alleen 


maar dichter naar evenwicht toe kan gaan, en niet in tegenovergestelde 


richting kan veranderen. De tegenstelling tot de Eerste Hoofdwet (die behoud 
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van energie uitdrukt) heeft de Tweede Hoofdwet een duidelijk statistisch 
karakter en is geldig voor het gemiddeld gedrag van een systeem. 

Een incidentele gebeurtenis, zoals een fluctuatie uit evenwicht hoeft 
niet precies aan deze wet te voldoen. Maar de kans op een ernstige schen= 
ding van de Tweede Hoofdwet is uiterst gering voor macroscopische syste 
men. Voor systemen bestaande uit slechts enkele deeltjes zijn de fluctu=- 
aties rond het gemiddelde gedrag relatief veel groter en de Tweede Hoofd=- 
wet verliest daar zijn betekenis. 

De entropie is een maat voor de wanorde van het systeem. Naarmate het 
systeem meer geordend is, is de entropie lager. En naarmate de entropie 
lager is, is de macrotoestand scherp gedefinieerd en kunnen we met meer 
zekerheid de microtoestand van het systeem aangeven. Ook kan men zeggen 
dat we meer informatie hebben over de microtoestand van het systeem, naar= 
mate de entropie lager is. In de evenwichtstoestand is de entropie maxi= 
maal, de informatie is minimaal, en we kunnen dan de toestand karakteri=- 
seren met een minimaal aantal grootheden E‚V‚N. Het begrip entropie speelt 
ook een belangrijke rol in de informatietheorie. 

We zullen zien dat Boltzmanns postulaten ons in staât stellen de tempera- 
tuur, de druk en andere thermodynamische grootheden te definieren in termen 
van de entropie, dus van de moleculaire structuur van het systeem. 

Dit hoofdstuk had betrekking op geisoleerde systemen. Als we in het 
vervolg spreken over de toestand van een macroscopisch systeem, dan be- 


doelen we de macrotoestand. 
IV. De Temperatuur 


Het begrip temperatuur hebben we nog niet gedefinieerd. In het vorige 
hoofdstuk kwam het hele begrip eigenlijk niet voor. In dit hoofdstuk zullen 
we het begrip temperatuur invoeren in termen van de entropie. We bekijken 
daartoe eerst een geisoleerd systeem, dat bestaat uit twee deelsystemen, 
gescheiden door een vaste wand die warmte doorlaat. De beide deelvolumes 
v, en V, zijn dus vast. De deeltjesaantallen in beide systemen zijn N, en 
N,. Het totale volume is V = Va + Vos en het totale deeltjesaantal 1s 
N= NN + N. De energieën E,‚ en B, van de beide deelsystemen zijn niet 


1 Ë 


constant, maar wel de totale energie E = E, + E. 


gehele systeem is gespecificeerd 


Een macrotoestand van het 


zodra we de macrotoestand van de 
- beide deelsystemen kennen en het=- 
zelfde geldt voor de microtoestand. 


Het totale statistisch gewicht @ 
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van een macrotoestand van het gehele systeem is dan hét product van de 


statistische gewichten Y, en oe van de macrotoestand van de deelsystemen. 


e= OA (h.1) 


Uit de definitie (3.11) van de entropie volgt dan 


ei (h.2) 


zodat de entropie additief is. 

Nu bekijken we de situatie dat elk van de deelsystemen in evenwicht is, 
maar dat ze nog niet in evenwicht met elkaar zijn. De toestand van Ast 
de deelsystemen is dan geheel bepaald door hun energieën E, en B, (immers 
Ns Vv, en Ns Vo zijn vast), waarvan de som constant is. De toestand van 
het systeem is dus bepaald door een enkele parameter, bv. E. 


De totale entropie als functie van E, 1S 


S= SEN + S(E - E‚) kid) 


o{ 
Er is evenwicht tussen beide deelsystemen wanneer de totale entropie S 


maximaal is. De evenwichtswaarde is dus 


a8 _ 
d8 0 (hl) 


1 
Volgens (h.3) wordt dit 
6 Ts) dS 
| +) = 0 (4,5) 


eni 
DE, 'V‚ oN, BE, ‘Vo oNo 


waar we expliciet aangeven dat we moeten differentiëren naar de energie, terwijl 
we volume en deeltjesaantal constante houden. Omdat E, =E= E, mogen we schrijven 
dS, 98 
(=S 


) ) 
BE, ‘VN, BE, ‘Vo oN, 


$ = 0 | 4.6) 
Als aan (h.6) voldaan is, zijn de beide systemen met elkaar in evenwicht, 
d.w.z. dat er geen warmtetransport meer plaats vindt van het ene deelsystgem 
naar het andere. 

We willen dat de temperatuur T een zodanige grootheid is dat geen netto 
warmtetransport meer plaats vindt tussen twee systemen met dezelfde temperatuur. 
Volgens (h.,6) bereiken we dat door de temperatuur te definiëren als een functie 


van ES. We definiëren nu 
dS. 


En = 


DE. von ST, 27 ê hik 
Led Ci 
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Ld 


De inverse temperatuur 1s gedefinieerd. als de afgeleide van de entropie 
naar de energie, bij constant volume en deeltjesaantal. De evenwichts= 
voorwaarde (h.6) zegt dan dat T, = T,. Met de definitie (h.7) is de tempe- 
ratuur identiek met de absolute temperatuur gedefinieerd met de ideale=gas- 
thermometer. Dat zullen we later neee: | 


Als nog geen evenwicht bereikt is, dan verandert de entropie volgens 


EE NE 
at at ° at 
E ® 
7 Eb, N nl iä EN N n 7 Gee le (4.8) 
EERE a Voolo 1 Zo 


Volgens de Tweede Hoofdwet moet de entropie toenemen, dus 


1 1 dE, | 


Als Ì, < Ís dan neemt E, toe, dus dan wordt er warmte getransporteerd van 


hogere temperatuur naar lagere temperatuur. Dit is in overeenstemming met wat 





we van het begrip temperatuur verwachten. 

De evenwichtsvoorwaarde (l.6) vonden we door te kijken naar twee deelsyste- 
men met vast volume en deeltjesaantal, die energie kunnen uitwisselen. 
We kunnen ook kijken naar twee deelsystemen met een beweegbare scheidings= 
wand. Dan zijn Vv, en Vo variabel met vaste som V = Vv, + V‚. Tenslotte kunnen 
we ook nog kijken naar de situatie van een scheidingswand die lekt, zodat 
de systemen deeltjes kunnen uitwisselen. (Dan kunnen ze overigens ook energie 


uitwisselen.) Als de systemen volume kunnen uitwisselen is er de evenwichts= 


voorwaarde 25 = 0 , waaruit volgt 
1 
6 Je; 
rh ps en N (4.10) 
|R: Be ee 2 2 | ’ 


Bij deeltjesuitwisseling vinden we de voorwaarde voor evenwicht 


d5, 95 

Gar) = () (h,11) 
oN, Es, AN, EV, | 

We zullen nog zien dat de druk van de deelsystemen in termen van de 


entropie is 
38. 


Bn 
av. 
L 


) is= 1,2 | (h.12) 


De 


JL 1 E. ‚N; ij 
zodat (h.10) samen met (h.5) betekentr- dat behalve de temperaturen ook de 
drukken gelijk moeten zijn in evenwicht. Ook kan men als in’ (h.9) inzien 


dat. het deelsysteem met grotere druk een toenemende volume heeft. 


| 
| 
r 
| \ 
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De chemische potentiaal u is gedefinieerd door 


es) 


= eh l = h.1 
1 1 5 


zodat H‚ = Ho de derde evenwichtsvoorwaarde is. Bij niet=evenwicht 


stromen de deeltjes van hoge naar lage chemische potentiaal. 


V. De Boltzmanmwerdeling. 


We willen Boltzmannspostulaten gebruiken om een waarschijnlijkheidsverde- 
ling af te leiden voor de microtoestanden van een systeem dat in thermisch 
evenwicht is in contact met een warmtebad met temperatuur T. Een warmte 
bad is een reservoir dat zo groot is dat de energie van het systeem veel 
kleiner is dan die van het warmtebad, zodat de temperatuur van het warmtebad 
niet beinvloed wordt door energieuitwisseling met het systeem. In evenwicht 
heeft het systeem dus een vaste temperatuur T, terwijl de energie niet con= 
stant is, maar kan fluctueren rond de evenwichtswaarde. De behouden groot- 
heden van het systeem zijn dus T, V,‚, N. (Voor een geisoleerd systeem waren dat 
E,‚ V‚, N.) De miecrotoestanden van het systeem nummeren we | is Ss aug En de 
bijbehorende energieniveaux zijn E‚s Es, «…… We zoeken de waarschijnlijk= 
heïidsverdeling over microtoestanden. 

De totale energie B van systeem plus warmtebad is behouden. Als het sys= 
teem zich bevindt in de microtoestand r met energie E‚, dan heeft het warmte- 
bad een energie E, = E‚. Het statistisch gewicht van deze energietoestand van. 
het warmtebad is Arp (Fo en FE). Omdat het systeem en het warmtebad samen een 
geisoleerd systeem vormen, geldt Boltzmanns eerste postulaat en de kans P, 
dat het systeem zich bevindt in de microtoestand r met energie E, 1S 


Py = CAE ij E) (Bel) 


De constante c is bepaald door normering 


DN P, = 1 (5.2) 


Op grond van de definitie van entropie (3.11) kunnen we (5.1) schrijven als 


Pp 


„=cexplS(E — E)/K) . (5.3) 


EE DE ME MENE NE OG ne TE 
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Omdat het warmtebad zoveel groter is dan het systeem, zijn de energieën E 


van het systeem die relevant zijn veel kleiner dan E . De waarde van E_ = E_ 


in (5.3) vertoont dan ook slechts relatief kleine fluctuaties rond een ge= 


middelde waarde En voor de energie van het warmtebad. We ontwikkelen de 


entropie in (5.3) in een Taylorreeks om Es 


35 (E_) 
— _ wb wb 
SE B) Sr Ero) 5 (E‚ B Ero) DE ij 
ln em en 3 Sol F) + 7 (5.4) 
E, O r wb es 5  s © Ld e 
AE 6 
wb 
Uit (h.7) weten we dat 
OE Ti 


waar T de temperatuur is van het warmtebad. De derde term in het rechterlid 
van (5.h) bevat dus de temperatuurverandering van het warmtebad tengevolge 
van warmteuitwisseling met het systeem. Omdat het warmtebad zo groot is, zijn 
deze en hogere=orde=termen verwaarloosbaar. We vinden dus voor de kans P, 

op de microtoestand r 


Sie Cpjole) (5.6) 


r 


waar @ = (Kr) Î, en 
Z(8) = sp (ST) 


De normeringsconstante Z wordt de toestandssom genoemd. De sommatie in (5.7) 
loopt over alle microtoestanden van het systeem. De temperatuur T is de 
enige eigenschap van het warmtebad die in de verdeling voorkomt. 

De belangrijke waarschijnlijkheidsverdeling (5.6) staat bekend als de 
Boltzmannverdeling. Zij is afgeleid voor een systeem In contact met een 
warmtebad.- We verwachten echter niet dat de toestand van een geisoleerd 
systeem zo geweldig veel zal afwijken van de toestand han hetzelf systeem 
in contact met een warmtebad. Bijvoorbeeld, een liter water in een calorimeter 
heeft dezelfde eigenschappen als een liter water in contact met de buiten= 
lucht van dezelfde temperatuur. De Boltzmannverdeling wordt dan ook veel | 
gebruikt voor geisoleerde systemen. 

De Boltzmannverdeling heeft vele plezierige eigenschappen. De verwachte 
waarde van de energie E van het systeem is de gemiddelde waarde van B 


Lr £ dB 


E= Zp E =-SdL | | A58 
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Om de entropie te vinden, hebben we een algemene formule nodig voor de 
entropie van een niet-geisoleerd systeem waarvan de toestand gespecificeerd is 


door de waarschijnlijkheidsverdeling P. Deze formule luidt 
S= „KEP, log p‚_. | | (5.9) 


We kunnen dit het gemakkelijkst afleiden door een groot aantal M identieke 
systemen te beschouwen, die alle verkeren in dezelfde toestand, gespecificeerd 
door de waarschijnlijkheden P,- Deze verzameling identieke systemen noemen 
we een ensemble. We beschouwen nu het statistisch gewicht van het ensembke, 

wanneer M, systemen verkeren in de microtoestand L n M, in 2, ete. Dit 
statistisch gewicht is weer gelijk aan het aantal manieren waarop het en= 
semble verdeeld kan worden in een groep van M, systemen, in een groep van 


M,, systemen, etc. Dus 


Q zite (510) 


Met Stirlingsformule vinden we dan voor de entropie van het ensemble 


S = k log Q 


= k [Mog M- ZM log M |} (5.11) 


aangezien Z M_ = M. Omdat de kans op een microtoestand r gegeven is door 
Ps is het aantal systemen in r juist p‚M. Fluctuaties rond deze waarde zijn 
mogelijk, maar als M erg groot is, zijn die verwaarloosbaar. De entropie 


van het ensemble is dus 


S= k [MilogM- ZE p,M log pM } (5.12) 
ofwel 
REN [M iogM-ZpMilogM=-Zp,M log p, 
En Rn Ë P, LOE Pi (5.13) 


Omdat de entropie additief is, vinden we uit (5.13) dat de entropie per 
systeem gegeven is door (59). Deze gelijkheid geldt voor een willekeurige 
waarschijrnlijkheidsverdeling. p 

In het speciale geval van de Boltzmannverdeling (5.6) is de entropie 


(5.9) gelijk aan 


e 


S = k log Ae LS 1d) 


a 


A TEE EE EE ECHEC vt Pe kak SiS mnd ek vraie: Sine 


EE Eg 


Oe ven ee de 
Ei 
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Tenslotte willen we het verband afleiden tussen de toegevoerde warmte 
en de entropie. Dat gaat het gemakkelijkst door eerst een kleine variatie in 
de gemiddelde energie E te beschouwen. Voor een willekeurige waarschijn= 


lijkheidsverdeling P. geldt 
dE = a(z PE) = EBdp, + ZE p,dE,. (5.15) 


De energleniveaux B hangen af van de vorm en de grootte van het volume V, 
maar niet van de waarschijnlijkheden P…: Dit is duidelijk als we ons herin= 
neren dat E‚_ de eïgenwaarden zijn van de totale moleculaire Hamiltoniaan. 


Dus bij een constant volume geldt dE, = 0, en dus 
dE = ZJ E dp. bij constante V . (5.16) 


Vanwege de definitie (h.7) van temperatuur mogen we schrijven voor evenwichts= 


toestanden 

ds = — dE + adV (5.17) 
(net deeltjesaantal N houden we constant) waarin a = G) . Met (5.15) 
volgt hieruit 

ds = = (2 p_dE, +2E dp) + adV (518) 


je 


Uit (5.9) is duidelijk dat dS = O wanneer de P, niet veranderen, dus wanneer 


dp, = OQ voor alle r. Dus 
En: | 
0 = FE p,dE, + adV, (5.19) 
en 
TdS = EE dp. (5,20) 


: In het speciale geval van de Boltzmannverdeling kunnen we (5.20) direct 
controleren, door te gebruiken 


L 
mn g Gog Pp, + log Z) (521) J 


De energieniveaux E‚ hangen alleen af van het volume, zodat we mogen schrijven 
dE 


Oh | (5.22) 


£ PoE, © 


Bij een reversibele adiabatische volutieverandering veranderen wel de E 


maar niet de p- Dus (5.22) geeft de gehele energieverandering bij een re= 


versibele volumeverandering. Volgens (2.6) is dE =-pdV bij een reversibele 


Del 


B a 


CT 
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adiabatische yolumeyerandering. Dus 
dE | 
== À ne 
à r Pr av (5.23) 


Uit (5.15), (5.20), (5.22) en (5.23) vinden we nu algemeen 
dE = TdS — pdV | (5.24) 


Dit is een algemene gelijkheid voor het energieverschil tussen twee dicht 
bij elkaar liggende evenwichtstoestanden van een systeem. Er komen alleen 
toestandsgrootheden in voor, zodat het geen verschil maakt of de overgang van 
de ene toestand naar de andere reversibel of irreversibel plaats vindt. 
Vergelijking (l.12) voor de druk volgt direct uit (5.24). 

In het speciale geval van de Boltzmannverdeling vinden we voor de druk 


p uit (5.23) en (5.21) 


_ 1 ò log 4 
Door (5.24) te vergelijken met (2.6) en (2.7) vinden we 
ds = oe voor een reversibele verandering , (5.26) 
as > 8 voor een irreversibele verandering . (5.27) 


T 


Dit geeft een direct verband tussen de entropie en de toegevoerde warmte. 
Uit (5.8), (5.1h) en (5.25) blijkt dat voor een systeem dat voldoet aan de 
Boltzmannverdeling de thermodynamische grootheden E‚ S, p geheel bepaald 
worden door log Z en de afgeleiden daarvan. Dit geldt heel algemeen. Het 


is gebruikelijk om in te voeren de Helmholtzvrije energie F door 


F=E- TS . | (5.28) 


Uit (5.14) volgt dan voor de Boltzmannverdeling 
1 
F == 5 log 2 (5.29) 
De relatie (5.24) kan dan ook geschreven worden als 
ar = - SaT - pdv , (5.30) 
waaruit-bijvoorbeeld volgt 
oF 


jen CE). ii (5.31) 


(vgl. (5.25)). Vergelijking (5.30)is praktisch voor een systeem in con- 


tact met een warmtebad, waar T en V de fundamentele variabelen zijn Lspes 


he ae en en: hanin 








Bij T= 0 (B 
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E en V. Uit (5.24) volgt 


== GE SR (5.32) 
Dus zoals E de effectieve potentiële energie is bij constante entropie 
(voor een geisoleerd systeem), zo is F de effectieve potentiële energie 
bij constante temperatuur (voor een systeem in een warmtebad). 

Het is niet ongebruikelijk om (5.26) te gebruiken als definitie van tem= 
peratuur, speciaal in puur macroscopische beschouwingen. Men zegt dan dat , 
er een grootheid T bestaat, zodanig dat voor reversibele processen ae ì 
de verandering is in een toestandsgrootheid. De temperatuur Is dan gede= 
finieerd als een "integrerende factor!" van de toegevoerde warmte &Qq. Wij 


geven de voorkeur aan de definitie (h.7), die minder abstract en meer aan= 


schouwelijk is. 
VI. Paramagnetische stof. 


We willen de resultaten van het vorige hoofdstuk toepassen op een een= 
voudig voorbeeld van paramagnetische kristallen. Dat zijn zouten bestaande 
uit grote moleculen (bv. Ce,Mg3 (NO), . 2 H_0), waarvan een deel een mag- 
netisch dipoolmoment heeft, tengevolge van het impulsmoment van de electro= 
nen. Doordat de moleculen zo groot zijn, liggen de magnetische dipool= 
momenten ver uit elkaar ei is de wisselwerking gering. Daarom beschouwen 

| we een modelsysteem, bestaande uit niet-wisselwerkende magnetische dipolen 
H‚ die in een uitwendig magneetveld B elk twee mogelijke energietoestanden 
hebben (dit komt overeen met impulsmoment 5), met energie =uB als de dipool 
parallel met het veld staat, en uB als de dipool antiparallel staat. We 
beschouwen het systeem in contact met een warmtebad met temperatuur T, zo= | 


dat we de Boltzmannverdeling mogen GAUSS Voor een enkele dipool is de 
toestandssom 


Z, es la + ene = 2 cosh (BuB) . (6.1) 


Omdat de dipolen niet wisselwerken, is een microtoestand van het hele systeem 


van N dipolen gegeven door de toestand van elk van de afzonderlijke dipolen. 


De toestandssom voor N dipolen is dan 


N 


Z(M,B,T) = Z, = [ 2 cosh (E45) ii | (6.2) 


ee 


De energie van hele systeem is 


9 log Z(N,B,T) 


FE En TR — NuB tanh (BUB). | | | (6,3) 


ĳ 


») is de energie gelijk aan -NuB. In de limiet van hoge tem- 
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peratuur (@ > 0) gaat de energie naar nul. De fracties p, en p_ van dipolen 


die parallel resp. antiparallel aan het veld staan bedragen 


B 


doma 
BuB „BHB 

À D= SD TKK (6.4) 
hee + 2 eosh (Bug ° *- 2 cosh (BuB) | 
f | Bij T = O0 is Pp, = Î, p_ = 0. Dus dan is het systeem volkomen geordend, en alle 
Rs dipolen staan parallel aan het veld. In de limiet van hoge temperatuur gaan P, 

en p_ naar 5, zodat dan de dipolen volkomen willekeurig geordend zijn. Een 
sterker magneetveld heeft precies hetzelfde effect als lagere temperatuur. 
| Het magneetveld, dat de dipolen probeert te richten wordt tegengewerkt door 
Je de temperatuurbeweging, die juist wanorde (= grotere entropie) probeert te 
hl veroorzaken. 
| Omdat dit model exact oplosbaar is, kunnen we uitdrukkingen vinden voor 
Ar alle fysische grootheden waarin we maar geinteresseerd zijn. De entropie S 
vinden we uit (5.1l), (6.2) en (6.3) 
ie | 
| S = Nk [ log (2cosh(8uB)) — BuB tanh (8uB) } (6.5) 
zn | De entropie is een monotoon stijgende funetie van de temperatuur T. Bij T = 0 
geldt S = O0, In de limiet van hoge temperatuur gaat S naar Nk log 2 = k log Ean 
| Dat is begrijpelijk, omdat bij T = O alleen de volledig geordende toestand 
| mogelijk is, terwijl bij T = ® alle toestanden mogelijk zijn. De warmte-capa- 
à citeit Co bij constant veld is de hoeveelheid warmte die nodig is om het 
a Ei | jj 
systeem bij constant veld een graad In temperatuur te doen stijgen. 
| Dus 
ik à DE 
ji es eq sz. __—_ n 
£ Cp Gra * Gora (6.6) 
Dn Immers, als geen arbeid wordt verricht, geldt dq = dE. 
$ Uit (6.6) en (6.3) vinden we 

2 1 
Cg = Nk (BUB) gn (6.7) 


sosh (BUB) 





De warmtecapaciteit is nul als T = O0, en gaat ook naar nul als T > @, 
Ongeveer bij de temperatuur waarvoor BuB = 1 heeft C5 een maximum. 

hs De overgang van volledige orde bij T = O naar volledige wanorde bij hoge 
temperatuur vindt feitelijk plaats over een eindig temperatuurgebied. 


Tenslòtte kunnen we het totale magnetische dipoolmoment M van het sys- 


| teem berekenen. We vinden 
__ 1 3 log Z(N.B.T ep 
n | M= Nulp, =P) =5 2108 ZINE) yy tanh(8u5). (6.8) 


WE VEREEN LOE EE en TN 
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Bij T = O is het totale dipoolmoment Nu. In de limiet van hoge temperatuur 


(of in de limiet van een zwak magneetveld) geldt 
2 
M= Nu BB | (6.9) 


zodat dan M lineair in B is. 


Uit (6.3), (6.5) en (6.8) kan men direct afleiden 


dE = TÄS — MAB … | (6.10) 


4 


Dit is het analogon van (5.24) voor dit systeem. De arbeid wordt nu verricht 
door B te veranderen. 

Voor hogere impulsmomenten 1, 3/2, 2, ... zijn 3, kh, 5 ....… orientaties 
van het dipoolmoment t.o.v. het veld mogelijk. Voor deze gevallen kan boven- 
staande theorie gemakkelijk. uitgebreid worden. 

Voor reêle paramagnetische zouten is het gemeten verband tussen M en B 
goed in overeenstemming met (6.8), of de generalisatie daarvan voor hogere 
impulsmomenten. 

De grootheden die we in dit hoofdstuk berekend hebben zijn alleen Het 
magnetisch deel van de energie, de entropie en de warmtecapaciteit. Ook de 
beweging van de deeltjes (roostertrillingen) dragen bij. Niettemin is een 
afzonderlijke behandeling van de magnetische bijdrage gerechtvaardigd, door- 


dat de koppeling van dipoolmoment en roostertrillingen gering is. 
VII. Het Klassieke ideale gas. 


We beschouwen een gas van niet-wisselwerkende moleculen. Dit noemt men 
een ideaalgas. Tevens nemen we aan dat de moleculen geen inwendige structuur 
hebben. De waarschijnlijkheidsverdeling over de verschillende energietoe- 


standen is gegeven door de Boltzmannverdeling. De kans op een energietoe- 


stand met energie B, is evenredig met exp(- BE). We willen de verdeling over 


de impulsen van de deeltjes weten. Daartoe moeten we weten hoeveel 
hd mogelijk zijn waarbij ben 1 een impuls heeft tussen D, en 

D. + BD deeltje 2 een impuls tussen Pe en jn + dba enz. Omdat de deeltjes 
nlet wisselwerken, is het voldoende om de verdeling van de energietoestanden 
over de impulsen voor één deeltje te vinden. Voor het gemak beschouwen we 
een rechthoekig vat met ribben Ls Ls L. Het volume is dan 


e= LLL, k _ EE et 
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Uit de quantummechanica weten we dat de x-component van de impuls P, 
zodanig moet zijn dat de De Broglie-golflengte van het deeltje een geheel 


aantal malen past in de afmeting L, van het vat in de x-richting. 
h 


waarbij n, een geheel (positief of negatief) getal is. Hetzelfde geldt 


voor p, en P‚- De toestandssom van een deeltje is dus 


2 2 2 2 
pn ee nn Bh n n n, 
benen nent A 
x Y ZZ em Lb, hb, L, | 


We gebruikten dat de energie van een impulstoestand met impuls D gelijk 


is aan Á- , waar m de deeltjesmassa is. We verwaarlozen nu het discreet 


zijn van de impulsen. Dan kunnen we de sommaties in (7.3) vervangen door 
integraties. Deze benadering betekent in feite dat we het gas niet quantum- 


mechanisch, maar klassiek behandelen. De benadering is geoorloofd mits 


2 2 


Eh ede of KT >>} 5 4 (7.h) 
emLx 2mLx 





Dezelfde voorwaarde moet gelden voor Ly en L,. De voorwaarde (7.h) eist dat 
de De Broglie-golflengte behorend bij een typisch thermische impulswaarde 
/ mT veel kleiner moet zijn dan de afmetingen van het vat. In reêle 
situaties is dat altijd het geval bij niet extreem lage temperatuur. 


In de klassieke limiet kunnen we (7.3) schrijven als een integratie over 


> R _X W . d 
Pp, met Pae E ‚ enz. We vinden 
v nT: de) onm. 3/2 
1,= 3 Spe on (145) 
h 2m Bh 


De laatste gelijkheid in (7.5) vinden we met behulp van de integraal 


| fax Ge (5 î | (7.6) 


Omdat de deeltjes onafhankelijk zijn, zouden we misschien verwachten dat 

de toestandssom Z(N,V,T) voor N deeltjes gelijk zou zijn aan bn juist 
zoals in (6,2). Maar we moeten rekening houden met het feit dat de deeltjes 
onderscheidbaar zijn, omdat ze identiek zijn. Dat betekent dat twee toe- 
standen die uit elkaar ontstaan door verwisseling van deeltjes, niet als 
verschillend, maar als een en dezelfde toestand beschouwd moeten worden. 
Door voor de totale toestandssom Ga te nemen hebben we een toestand 

> > 


n 8 > > > > 
met impulsen en coordinaten (ps Pos «« Ds Fjs Fos «° En) net zoveel 








2d, 
malen geteld als er manieren zijn om N deeltjes onderling te verwisselen, 
Dit aantal is N!, en we vinden voor de toestandssom van een klassiek 


ideaalgas van N deeltjes 


af r2mm, 3N/2 
2 
\ Bh 


Z(N,V.T) = (zo) = | (7.7) 


Ke Ni 


| 


De factor (wi) kwam niet voor in het model van N magnetische dipolen in 
(6.2), omdat daar de deeltjes op vaste plaatsen in de ruimte zitten, en der- 
halve twee toestanden die uit elkaar ontstaan door zowel de orientatie als 
de positie van twee dipolen te verwisselen niet afzonderlijk geteld zijn 
in de toestandssom. 

De verdeling over toestanden kunnen we in de klassieke limiet schrijven 
als een continue waarschijnlijkheidsverdelingsfunctie f Er ee D) over 


impulsen, De functie f is gedefinieerd door de eis dat f Dn Pere DpldD, 


äp, de kans is dat deeltje 1 een ee heeft tussen D, en D. + dp,» 


deeltje 2 een impuls tussen B en Po je dp, enz. Uit de (genormeerde) 


Boltzmannverdeling volgt dan dat 


> > 2 


u HD zech) JT.) 


exp [ - 


PD, - Dy) = (oen 


De verdelingsfunctie van één deeltje is dan 


mn 
2D) = GE P ep (EE) (7.9) 


Deze verdeling staat bekend als de Maxwell-verdeling. Zij kan geschreven wor- 
den als het product van een verdeling over p‚ ‚ een over p en een over Pp, . 
Dit toont aan dat de componenten van de Impuls onafhankelijk verdeeld zijn: 
als je weet wat de x-waarde van de impuls is, dan zegt Je dat niets over 

de te verwachten y-waarde. De verdeling van elke component heeft de bekende 
klokvorm van de Gaussische verdeling. 

Een veel gebruikte verdeling is die van de absolute waarde van de impuls 
van een deeltje. De grootheid g(lp)dp is gedefinieerd als de waarschijnlijk- 
heid dat de absolute waarde van de impuls p = / + Je + p‚° een waarde 
heeft tussen p en p + dp. Hieruit vinden we direct 

2 
sp) = bep GEIN tE (7-10) 
omdat het volume van een bolschil met binnenstraal p en buitenstraal p + dp 


gelijk is aan rp“ As … ú 


Uit de toestandssom (7.7) kunnen we zoals gewoonlijk alle thermody- 


namische grootheden af leiden. 


De Energie E is 


Verre TE 


me 


aka 


bk kad 
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E= - SEE 2 yer (7.11) 


Voor een klassiek ideaal gas hangt de energie dus niet af van het volume, 
alleen van de temperatuur. Elk van de 3N bewegingsvrijheidsgraden draagt 
dus à kT bij tot de energie, | 

De druk p is te vinden uit (5.25) 


PV = NKT (7.12) 


Dit is precies de ideale-gas-wet, waaruit blijkt dat de definitie (h,7) 
van de temperatuur precies overeenkomt met de absolute temperatuur volgens 
de gasthermometerschaal. 


De warmtecapaciteit blj constant volume ce is gedefinieerd als de hoe- 


V 
veelheid warmte die nodig is om bij constant volume het systeem 1 graad in 


temperatuur te doen stijgen 


_ (3 _ 3E / 
Uib (Tell) rotet 
C = 3 Nx (Tk) 
Vv 2 3 d 


De warmtecapaciteit is dus onafhankelijk van volume en temperatuur en on- 
afhankelijk van het soort gas. | 
De entropie is te vinden uit (5.,1h), waar we voor log N! Stirlings 


formule gebruiken. Dan vinden we 


S = Nk [2 log( SE) + 3 + log D) ) A (7.15) 
Bh ‘ 
De entropie is dus evenredig met het deeltjesaantal N, en hangt op een= 
voudige wijze af van de dichtheid en massa. 

Als toepassing van deze uitdrukkingen voor E‚ p, en S en als illustratie 
van het verband tussen verrichte arbeid, toegevoerde warmte en entropie= 
verandering bekijken we nu de compressie van een ideaalgas. We bekijken het 
geval dat die compressie adiabatisch (&Q = 0) of isotherm (AT = O0) geschiedt. 
Verder onderscheiden we oneindig langzame (reversibele) en snelle (irre= 


versibele) compressie. Het beginvolume is V 


2 1 


1” het eindvolume is V‚ < V.. 


A. Isotherm en reversibel. 


De verrichte arbeid bedraagt nu 


Vo 
MW== {PdV , | OO CT. 16) « 
| 
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enn me menemen eer ee ek renee me 


EEEN MERE AT CORET EON MERODDAEED LOEDER OOKE MAESEN it VE TO TOER EE SON ME GAMOG TE OGER EEP ENNE 


RER OEE LN ee LE PE ren CRN ME ELNNEE ST Aint OE MT Arden en ee oe meh edn een ern en an edere RA nd es 





2, 


Ld 


wegens de reversibiliteit. De temperatuur T is constant, en pV = NKT. 


Dus 
Vv Vv | 
AW= = {TEE av = NKT log pe > 0 (7.17) 
1 2 | 


Dus er is arbeid op het systeem verricht. Doordat de temperatuur constant is, 


is dé energie wegens (7.11) onveranderd. 


AE =O (7.18) 


Uit de Eerste Hoofdwet volgt dan 


Vv 4 A 
AQ = = AW = = NKT log Ee, < 0 ' (7.19) 
2 


Dus het systeem heeft warmte afgestaan aan het warmtebad. Omdat het proces 
reversibel is volgt de entropieverandering uit (5.26) 


Vv 
AS = = Nk log (Ei) <0 (7.20) 
2 


De entropieverandering is ook te vinden uit (7.15). 


De entropie is dus afgenomen. 


B. Adiabatisch en reversibel. 


Nu geldt AQ = O0, zodat volgens de Eerste Hoofdwet 


Vo | 
AE = AW= = f pdV (7.21) 
1 


De temperatuur is nu echter niet constant. Volgens (7.11) gelät 


dE = 2 Nk dT (7.22) 


Op elk moment is het systeem in evenwicht, zodat 


AE = NAT = - pav =- EE av (7.23) 
Dus 
d(log T) = - Za (tog Vv) (7.2k) 
en 
T V | | 
2 1,2 
log (==) = log (—) /3 (7.25) 
EL Vv e 
1 | e 
waar T, de begintemperatuur en T, de eindtemperatuur is. Uit (7.11) volgt 
nu | 
E, T 


4 T 


V 
ES _ 2 (—)2/3 
1 0 


(7.26) 
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Ad 


Dus zowel de energie als de temperatuur stijgen bij compressie. De ver= 


richte arbeid AW is gelijk aan het energieverschil 


V 
jn pe nn: í ES Je NN 
M=E Ez LG) lk (7.27) 
Uit (5.26) volgt 
AS= 0 (7.28) 
Ook met behulp van (7.15) kan dit gevonden worden (ga dit na). 


C. Isotherm en irreversibel. 


Nu comprimeren we het gas met een eindige snelheid. Volgens (2.5) vol= 
doet de verrichte arbeid aan de ongelijkheid 


Vo LEEN 
AW > = { pdV = NKT log (o) (7.29) 
1 2 


Dus nu wordt meer arbeid verricht dan in geval A. De energie is onveranderd 
omdat T constant gehouden wordt. Volgens de Eerste Hoofdwet voldoet dan 
de toegevoerde warmte aan de ongelijkheid 
V | 
1 
AQ < = NKT log (—”) (7.30) 


Vo 


Nu wordt dus meer warmte afgevoerd naar het warmtebad dan in geval A. De 
begin= en eindtoestand zijn hetzelfde als in A, dus ook geldt 
V, 
AS = = Nk log (—) … CT) 


hi 


Deze entropievermindering voldoet aan de ongelijkheid (5.27). 


D. Adiabatisch en irreversibel. 


De verrichte arbeid is nu groter dan in geval B. Omdat AQ = O, moet 
AE groter zijn dan in B. Dus ook de temperatuur wordt meer verhoogd dan 
in B het geval was. De entropie=verandering AS is nu positief, in overeen= 


stemming met (5.27). 
VIII. Het klassieke reêle gas. 


Het ideale-gas-model is een benadering die alleen goed is bij lage dichtheid. 
In het kenden moeten we echter de wisselwerking tussen de deeltjes in | 
rekening brengen. De toestandsvergelijking die het verband geeft tussen dicht= 
heid, druk en temperatuur heeft niet meer een algemene eenvoudige ge= | 
daante, die universeel is voor alle Gsen zoals in (7.12) het geval was. 
Bij niet al te hoge dichtheid is een reeksontwikkeling van thermodyna= 


mische grootheden in machten van de dichtheid n = N/V nuttig. Zo'n reeks= 
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ontwikkeling wordt een viriaalreeks genoemd. De viriaalreeks voor de druk 


kan geschreven worden als 
p=nkT (1 +nB(T) +nc(T)+.….] (81) 


De ideale=gas=wet (7.12) is de eerste=orde-term van (8.1). 

Op grond van algemene kwalitatieve overwegingen kan men ook wel trachten 
een verbeterde toestandsvergelijking te gissen. Met groot succes is dit in 
de vorige eeuw gedaan door J.D. van der Waals. Hij ging uit van het idee 
dat de intermoleculaire wisselwerking een aantrekkend lange-drachtsstuk 
en een afstotend korte-drachts-stuk heeft. De harde afstotende kern van de 
wisselwerking verhindert dat de deeltjes elkaar overlappen. Daardoor hebben 
de deeltjes niet het totale volume V tot hun beschikking. Per deeltje moet 
van het totale volume een volume b worden afgetrokken, waarin b het effec= 
tieve volume van een deeltje voorstelt. Dus in de ideale=-gas=wet vervangen we 
V door V = Nb, Het aantrekkende stuk houdt de deeltjes enigszins bijeen. 

Dus er is meer energie nodig voor een deeltje om het gas te verlaten dan in 
een ideaal gas het geval zou zijn. Anders gezegd, de druk op de wand is niet 
zo hoog als voor een ideaal gas het geval zou zijn. Dat betekent dat we in 
de ideale=gas=wet de druk p moeten vervangen door de druk plus een positieve 
term. Deze extra term moet kwadratisch zijn in de dichtheid, omdat het to= 
tale aantal wisselwerkingen evenredig is met het kwadraat van het aantal 
deeltjes. De resulterende vergelijking is de toestandsvergelijking van 


Van der Waals 
2 
(p +n a)(1 - nb) = nkT (8,2) 


In de afgelopen vijftien jaar zijn modellen ontworpen, waarin (8.2) wordt 
gevonden als resultaat van een microscopische theorie. 

Het verband tussen p en V bij vaste T is volgens (8.2) aanzienlijk ge- 
compliceerder dan volgens de ideale=gas=wet. Als functie van V blijkt p 
lokale uiterste waarden te hebben, tenminste zolang T beneden een zekere waarde 
blijft. Dit wordt geïnterpreteerd als condensatie. Als we een gas iso- 
therm samenpersen, zal de druk toenemen naarmate het volume kleiner wordt, 
totdat de verzadigingsdampspanning bij de betreffende temperatuur bereikt 
is. Als we het volume nog kleiner maken, zal de druk eerst niet verder toe= 
nemen, maâr het gas condenseert tot vloeistof. Pas wanneer het volume zo | 
klein is geworden dat alleen vloeistof aanwezig is, zal de druk verder toe= 
nemen. In het pV-diagram zal de re ans een horizontaal stuk hebben. 
Condensatie vindt slechts plaats wanneer de temperatuur beneden een kri= 


tische waarde Ln ligt. 








af. 


Meri- zegt nu dat de beide lokale uiterste waarden in de pV-curve van 

Van der Waals geïnterpreteerd kunnen worden als condensatie. De golving 

in de pV-curve zou dan vervangen moeten worden door een horizontale 
lijn op zodanige 
hoogte dat de 
beide resulterende 
oppervlakjes ge- 

* Lak zt Det 

| horizontale stuk’ 

stelt dan con= 

Lo 

densatie voor. 

De golving wordt 

dus niet als reëel 


os beschouwd (tenzij 





in termen van me= 
tastabiele toe= 
standen, bv. on= 
0, to ‚Ss v* derkoelde damp). 
| De kritische tem= 
peratuur In volgens 
Van der Waals is dan 
de temperatuur waar 
de pV=curve nog Juist 
een stationair punt vertoont. De bijbehorende druk Pp, en dichtheid n, worr 
den de kritische druk en de kritische dichtheid genoemd. Uit de Van der 


Waalsvergelijking (8.2) kunnen P.» nh, En ha gevonden worden m.b.v. de eis 


Ë 


dat p als functie van V (of n) een buigpunt vertoont 


e 
SD ò Pp — | 
ov | 
Hieruit volgt: 
„8e. en A EN 
in 7 27 kb ° Ze  3b ’ Po 7 27b2 (Bu) 


De Van der Waalsvergelijking (8.2) kan in dimensieloze vorm geschreven wor- 


den 


Eea esn vaan 


3 (8.5) 


waar 
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p= P/P 3 n= n/n, 3 Tr = B Li (8.6) 


De figuur laat de Van der Waalscurven in dimensieloze vorm zien. Hier is 
* | * “ | 
V == V/V, en Le = N/n zodat V = 1/n … 
De overeenstemming met experimentele resultaten is redelijk bij niet te 


hoge dichtheid. 
Een meer algemene, moleculaire beschrijving van een niet-ideaal gas vin- 


den we door uit te gaan van een moleculaire Hamiltoniaan, bestaande uit 


de kinetische energie van de deeltjes, en een potentiële energie . 
vr, Ee RE die een functie is van de positie van alle N deeltjes. De 


eenvoudigste aanname is dat V geschreven kan worden als de som van wissel= 
werkingspotentialen voor paren van deeltjes, waarbij de paarwisselwerking 


alleen afhangt van de afstand tussen de twee deeltjes 


en > > > > 
V(r,, Tse Ty) = 135 v(| ri =r; Ek 2 (8.7) 


De verdelingsfunctie over de impulsen en de posities van alle deeltjes 
Co Danae Ps Es Bas Ea) kunnen we vinden met een soortgelijk argu= 


ment als waarmee we in hoofdstuk III de Maxwell=Boltzmannverdeling vonden 
voor een ideaal gas in een uitwendig veld. Men kan dit zelf nagaan door 


een ensemble van identieke systemen te beschouwen, en het statistisch ge= 


wicht voor een verdeling van microtoestanden Dn Es 
>> 


r.)te 
| n ee . e N 
beschouwen. De verdelingsfunctie Ê(Dys---Dy5 PERREN is dan 


> za. > 5 1 4 =B(K + Vv) 
FCDys «Ds Eys <-Fy) = Z(N,V,T) N!n3N © ‚ (8.8) 
waar K de kinetische energie is 
a 
bj 
K = L Dm (8.9) 


en waar V de potentiële energie is, gegeven door (8.7). De toestandssom 


Z(N,V,T) is een generalisatie van het ideale=gas=resultaat (7.7) 


zE) 





ANN en GED äpy är, rin (8.10) 
| N'h 


Een typische paarwisselwerking vr) heeft een afstotend kortesaraehtestuk 
en een aantrekkende lange=drachtsstuk in de figuur. | 
Vergelijking (8.10) voor de 
toestandssom dient als uitgangs= 
punt om thermodynamische groot= 
heden te vinden, bijvoorbeeld 
de druk p(N,V,T). Een algemene 
eenvoudige gedaante voor de 


toestandsvergelijking 
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is niet meer te vinden en we moeten benaderingen maken. Wel is het mogelijk 
om expliciete uitdrukkingen te vinden voor de viriaalcoefficienten B, C, 
ete. We zullen dit laten zien voor B(T) (de tweede viriaalcoefficient). 
Omdat de Hamiltoniaan de som is van de kinetische energie (die alleen van 
de impulsen afhangt) en een potentiële energie (die alleen van de posi= 
ties afhangt), kunnen we de toestandssom (8.10) schrijven als het product 


van de ideale-gas-toestandssom (7.7) en een correctiefactor, 
Z(N,V,T) = Zig Wp > (8,11) 


waar Gr de configuratie-integraal is 


Ar = se if A IN A | (8.12) 
vi 1 N | 
ienmege: (Ba): Bnei ne Ld schrijven als een product van paarwisselwer= 

kingsfactoren 
hes > 
-BV _ =Bv(|r. - r.|) 
e = 3; e 3 J R (8.13) 
We voeren nu in de functies his als 
> > 
1} 
> 


Dan gaat EE naar nul als (r, _ r,) > oo, We kunnen On schrijven als 


Qq en AK} > > | 
N S dr. dr, Il + Az 3) 


ln % 
; 


mn 
=! Arg: arg [1 + D,s + EE), ; Ag + er J (8.15) 
=1 + rl far Alr) + .… 
De Helmholtz vrije energie F is dan 
5 Ee 4 
F= F3 6 log Wp - (8.16) 
Voor de druk vinden we dan uit (5.25) 
_ U is 
P = Pia* B âv [: MNI) fer 1) + ee (8.17) 
„NKT -— ZKT N2 f > N3 | 
5 D Jer) + eD | (8.18) 


Men kan direct nagaan dat de verdere, termen van hogere orde in de dichtheid 


BT Me 


N 
À 
Á 
8 
| 
, 
b 
f 
8 
Ì 
| 
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„De term 1-e 


Door dit te vergelijken met (8,1) vinden we 


B(T) = If Mr) = on jr de [rev (8.19) 


-Bvlr) is ongeveer 1 als v(r) groot is (voor kleine r), 


en is ongeveer Bv(r) als v(r) klein is (voor grote r). We benaderen nu 


j-eBv(r) = 1 als o<r Fo (8,20) 


R 


Bv(r) als r>r- (8.21) 


waar Fo de afstand is waar v(r) = o. Dan is Yogeen maat voor de diameter 


van de deeltjes. Door (8,20) en (8,21) in te vullen in (8.19) vinden we 


_B(T) 5 r + 2rB | dr r vr) . (8.22) 
5 | 


In deze simpele benadering bestaat de tweede viriaal-coefficient uit twee 


termen, êên die onafhankelijk is van de temperatuur, en één evenredig met 


E ‚ Dit kunnen we vergelijken met de viriaalontwikkeling van de Van der 


Waalstoestandsvergelijking (8.2) 


p= at Úi+n(b-e) + + … ] | (8.23) 


Door verlijking van (8,22) en (8.23) vinden we 


B (8.2) 
3 ©) 
a = = Ôm Á är ee ete | (8,25) 
/ | 


© 


Dit geeft dus een moleculaire interpretatie van de grootheden a en b. 
Inderdaad heeft b te maken met het eigen volume van de deeltjes, en a 


met de lange-drachts aantrekkende wisselwerking. 


Dit is een voorbeeld van een eenvoudige berekening van meetbare grootheden 
(B(T), a, b), in termen van de wisselwerkingspotentiaal v. Ook voor meer 


gecompliceerde beschouwingen is (8,10) een geschikt uitgangspunt. 











IX. Systemen met deeltjesuitwisseling. 


In hoofdstuk III beschouwden we een volledig geïsoleerd systeem, waarvoor 
E, V, N vaste waarden hebben. Daarbij behoorde de verdeling over micro- 
toestanden, waarbij alle microtoestanden met deze waarde E‚V‚N een even grote 
waarschijnlijkheid hadden. Deze verdeling wordt wel de microcanonieke verde- 


ling genoemd, In hoofdstuk V beschouwden we een systeem in contact met een” 


_warmtebad (dus met fluctuerende energie) maar met vast deeltjesaantal MN 


Hier hebben T,‚V‚N vaste gegeven waarden. De bijbehorende verdeling is de 
Boltzmannverdeling (5.6), die ook wel de canonieke verdeling genoemd wordt. 

In dit hoofdstuk beschouwen we een systeem dat niet alleen energie, maar 
ook deeltjes kan uitwisselen met een groot systeem. Het warmtebad fungeert 
dan ook als deeltjesreservoir. Het systeem kan bijvoorbeeld omgeven zijn door 
vaste, lekkende wanden, of het kan een denkbeeldig hokje zijn in een groot 
systeem, bijvoorbeeld een kubieke centimeter in een gascylinder. Dan kunnen 
we vragen naar de waarschijnlijkheid Py ‚r° dat het systeem N deeltjes bevat, 
en dat dit N-deeltjessysteem verkeert in de microtoestand r. Het systeem met 
het warmtebad is in evenwicht. De afleiding gaat hetzelfde als voor de Boltz- 
mannverdeling (zie hoofdstuk V) alleen met dit verschil dat niet alleen 
de energie E_, van het systeem maar ook het deeltjesaantal N een variabele is, 
De energie E‚ en het deeltjesaantal N, van systeem plus warmtebad is constant. 
De entropie van het warmtebad Sp Eo En NN ) kunnen we ontwikkelen in een 
dubbele reeks. Het resultaat is dat 

Py,r cPEnr SEN ne(r,von) (9.1) 

waar B = (kr)”!, en waar T de temperatuur en u de chemische potentical van het 
warmtebad is. We gebruikten de definitie (h.13) van chemische potentiaal. De 
verdeling (9.1) wordt wel de Gibbs-verdeling of de groot-canonieke verdeling 
genoemd. De normerings-constante E(T,V,u) wordt de groot-canonieke toestandssom 
genoemd, terwijl de "gewone!" toestandssom Z(T,V‚N,), (5.7) ook wel de canonieke 
toestandssom wordt genoemd. De groot-canonieke toestandssom Z(T,Vv‚u) 1S gedefi- 
nieerd als 


Z(T‚v‚) = È Ee PEnrtBeN _Ì BEN z(r,v‚n) … (9.2) 
N=o r „ N=o | | 


32. 


De Gibbs-verdeling is gespecificeerd door de drie vaste waarden T, V, u. 


Voor vast deeltjesaantal vonden we (5.2l) 


1 B 
= — + 
tn (9.3) 
voor het verschil van twee evenwichtstoestanden. 
Als ook N kan variëren, dan wordt dit 
| P } | 
AS =GdE+EV- rd (9.h) 
waar we (h.13) gebruikten. Dus bij variabel deeltjesaantal wordt 
(5.2) gewijzigd in 


dE = TdS — pdV + uaN (9.5) 


6 E 
Verg. (9.4) of (9.5) geeft een relatie tussen veranderingen in thermo- 
dynamische grootheden, voor twee dicht bij elkaar liggende evenwichts- 
toestanden, 


Voor de Gibbs-verdeling geldt volgens (5.9) 


D = =kl Pr. log p 
Ne Nr Nr (9.6) 


= _k] BuN BE — log % J 
NrPNr le 


—uN E 
Ered T + k log X . 


waar N en E het gemiddelde deeltjesaantal en de gemiddelde energie zijn. 


We definiëren nu de zogenoemde grand potentiaal Y door 


y= 5 log X (T, V‚u) | (9.7) 


Uit (9.6) volgt dan 


V=E-TS -uN= F- UN | (9.8) 


waar F de Helmholtz vrije energie is, gedefinieerd in (5.28). De grand 
potentiaal Y speelt in het geval van de Gibbsverdeling dezelfde rol als 

de Helmholtz vrije energie in het gewal van de Boltzmannverdeling, zo- 

als blijkt door (5.29) en (9.7) te vergelijken. De overige thermodynamische 
grootheden volgen nu door differentiatie van log 4(T, menen de 


krakteristieke variabelen 





* 


33. 


En Nn do Z(T Vk ) ks En 9 \ 
B =P Pr Ene Cg GAT ie ET | zE 9 
: ks ZV 
e _ 1 (2log Z(T,V,E), | j 
r 
_ À (log Z (T,V,P), | 


T‚u 





De gelijkheden (9.9) en (9.10) volgen rechtstreeks uit (9.1), 
(9,11) volgt het gemakkelijkst door (9.8) en (9.5) te combineren 


AY = dE — TâsS- SAF — uäN — Nau 


=-pdV — SäT -— Näu „ | | (9,12) 


waarmee (9.11) triviaal is wegens (9.7). We kunnen ook schrijven 


ee) 


Pp == ls | 
V (9,13) 


waaruit blijkt dat de grand potentiaal de effectieve potentiele energie 
is bij constante temperatuur en chemische potentiaal (dus voor een systeem 
dat energie en deeltjes kan uitwisselen met een warmtebad). 

Door vergelijking zal duidelijk zijn dat de thermodynamische inter- 
pretatie van de microcanonieke verdeling, de canonieke verdeling (of Boltz- 
mannverdeling) en de groot-canonieke verdeling (of Gibbs-verdeling) 
eenzelfde structuur heeft. De relaties tussen thermodynamische grootheden 
is identiek in de drie gevallen, maar in elk van de gevallen hebben we een 
andere drietal karakteristieke grootheden. De microcanonieke verdeling heeft 
als karakteristieke grootheden E‚V‚N. De relatie van destatistische verdeling 


met de thermodynamica is de relatie (3,11), 


S = k log Q | | (9.1L) 


die de entropie S geeft als functie van E‚V‚N, Andere thermodynamische 
grootheden volgen door differentiatie van S naar deze karakteristieke 


grootheden, zoals blijkt uit de algemene thermodynamische differentiaal (9.4) 


8 | 
(oe) os os 

„= dE ‚, P=T (5 ue T (5e) 

T N,V 9V NE oN E‚v (9.15) 


e 








3h , 


De canonieke (of Boltzmannverdeling) heeft als karakteristieke 


grootheden T,V‚N. De relatie met de thermodynamica is (5.29) 
F = 5 log Z , (9.16) 


die de Helmholtz vrije energie geeft als functie van T,V‚N. De 
differentiaal dF is uitgedrukt in de veranderingen in de karakteris- 


tieke grootheden door 
aF = -SäT -— pdV + uäN , | __ (9.17) 


zoals men vindt uit (9.5) en 65,28). Andere thermodynamische groot- 
heden vinden we door F te differentieren naar de karakteristieke groot- 
heden 


en ren sey «HO 


Voor de groot=canonieke verdeling (of Gibbs-verdeling) zijn de 
karakteristieke grootheden T, u, V. De relatie met de thermodynamica is 


(9.7) 
W = = log X (T,V,u) | | | (9.19) 


Andere thermodynamische grootheden vinden we door (9.19) te 


differentieren naar de karakteristieke grootheden. Uit (9.12) volgt 


gy 9y | gy | | 
p=- (5) Ee Bea Go N= e= Gr n (9,20) 


We willen er de nadruk op leggen dat de relaties tussen thermodynamische 


grootheden (zoals (9.h), (9.5), (9.8), (9.12), (9.13), (9.15), (9.17), (9.18), 


(9,20)) algemeen geldig zijn voor systemen in evenwicht. | 

In het bijzonder zijn ze alle geldig voor de drie evenwichtsverdelingen 
die we besproken hebben. De drie verdelingen gelden echter voor verschil- 
lende situaties, en de meest natuurlijke relatie met de thermodynamica is 


verschillend voor de drie verdelingen. 


e 
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|. De warmtecapaciteit van een gas bij constant volume Cy Is gedefinieerd als 


- (Q 
vG | 
De warmtecapaciteit bij constante druk is 
dq 
c_ = (&) 
AT n 
i Pp 
waarbij de (eventuele) volumeverandering reversibel is. 4 ’ 
a. Bewijs dat 
C. = SE 
V AT 
dE 
pT GP + py 
Pp ® 


waar yde thermische uitzettingscoefficient EE IS. 
P 


2. Voor een ideaal gas geldt 
pV = NkT 


en 


me: 


a. Bereken cy en ep voor een ideaal gas. 


b, Laat zien dat bij een adiabatische reversibele toestands- 
verandering van een ideaal gas geldt 


C‚dT = -pdV 


V 
Laat zien dat voor. zo'n verandering geldt 
C, e =C | 
T N vP Ne constant 
ee | | 
p(v P 5 = constant 


3. De enthalpie H van een systeem is gedefinieerd als 





H=E +pV . Druk Es uit in H. 











een gien: ne ee 


4, 


We beschouwen het volgende eenvoudige model van een kristal van N atomen. 
Elk atoom kan zitten op een roosterpunt met energie nul, ofwel op een 
plaats midden in de roostercel met energie e. Er zijn dus N roosterpunten, 
en N alternatieve roosterpunten. Als een atoom niet op zijn roosterpunt 
zit, kan het op elk van de N alternatieve posities zitten. 

De atomen zijn niet onderscheidbaar. Beschouw nu een toestand waarin n van 
de N atomen niet op hun roosterpunt zitten. 


ES 


a. Wat is het statistisch gewicht van deze toestand? 


b. Maak de Stirlingbenadering voor de faculteit log p: = p log p = p. 
Geef de entropie en de temperatuur van deze toestand. 


c. Hoe groot is de factie 5 als e = leV, T = 300K? 


| 
Zo eV als T = 300 K). 

Beschouw een ideaal gas in een geisoleerd bolvormig vat. Bij botsingen met 
de wand verandert het impulsmoment t.o.v. het middelpunt niet, zodat het | 
totale impulsmoment een behouden grootheid is. Beschouw nu een toestand met 
energie E en totaal impulsmoment L . 


1 (KT = 


a. Geen een uitdrukking voor de verdelingsfunctie £(p‚t) van de deeltjes, als 
functie van E en L . 


mrs 
b. Wat is deze verdelingsfunctie als L = o? 


De warmtecapaciteit van 1 gr. water is 4,2 J/K. 

We brengen Ì kg water met een temperatuur van 20 C in contact met een 
warmtebad van 80 C. We wachten totdat zich het evenwicht heeft ingesteld, 
en het water ook op 80°C is gekomen. 


a. Wat is de entropietoename van het water plus het warmtebad? 
Nu verwarmen we het water in twee stappen, door het eerst in contact te 
brengen met een warmtebad van 50 C, en vervolgens met een warmtebad van 
80 CG: ’ 


b. Wat is nu de totale entropietoename? 


c. Wat wordt de entropietoename als we het aantal stapjes naar oneindig 
laten gaan? 


J. We bekijken een systeem van N harmonische oscillatoren, elk met energie 
2 eigentoestanden bij de energieën E‚ =he® (n + }). 


De oscillatoren zitten vast in de ine 


a. Geef de Helmholtz vrije energie F, de gemiddelde energie E, 





de warmtecapaciteit c, de entropie S en de chemische potientiaal u 


Schets deze grootheden als functie van de temperatuur en als functie van W. 


Dit is het Einstein=model voor de roostertrillingen van een kristal. 


gg 


8. We beschouwen een tweetal cylinders met gas, elk afgesloten door een 
zuiger. De zuigers zijn gekoppeld viaeen hefboom als in de tekening. 
De oppervlakte van de ene zuiger is O,, van de andere 0. | 
De armen van de hefboom hebben lengte a,, en a 
De systemen kunnen tevens warmte uitwisselen. 


a. Geef de voorwaarden voor evenwicht. - 





draaipunt 





9. Een geisoleerd systeem is door een beweeglijk adiabatische wand in twee 
deelsystemen verdeeld. De hl kunnen geen warmte uitwisselen, maar 
bij een volumeverandering wordt onvermijde- 
lijk ook energie uitgewisseld. Dus volume- 
verandering en energieerandering.zijn nu niet 
2 onafhankelijk. Immers 





dE; 2 -p‚dV, | (1) 
Volgens de fundamentele thermodynamische differentiaal geldt 


dE, = Tds, — p;dV, (ii) 





ĳ 
| 
ü 
| 
| 
E 
| 


Verder zijn de totale energie E‚ + E, en het totale volume constant. 
De temperaturen zijn verschillend. | 


a. Laat zien dat bij een spontane verandering geldt 
(pj =P) dV, Zo , 
zodat het volume toeneemt waar de druk het grootste is. 


b. Laat zien dat dS = o, zodra Pj * Pj- Dus bij evenwicht zijn de drukken 
gelijk, ook al zijn de temperaturen verschillend, 


- 4 — 


10. Voor een paramagnetische stof van spins } is volgens (6.5) de entropie 
een functie van BuB. | 
Neem aan dat we het magnetisch veld B adiabatisch en reversibel reduceren 
tot de lagere waarde B'. | 


a, Wat is de eindtemperatuur T' ? 
(Deze methode van adiabatische demagnetisatie wordt gebruikt om 
materie dicht bij het absolute nulpunt nog verder af te koelen.) 





| Il. Een ideaal klassiek gas van N deeltjes is opgesloten in een deelvolume Vv, 
van een groter volume V. 

ĳ Het gas wordt door een wand in V, gehouden. 

Op zeker moment wordt deze wand plotseling weggetrokken, zodat het gas 


het gehele volume V gaat vullen. Dit is een irreversibel proces. 


| a. Bereken de verandering in energie, alsmede de toegevoerde warmte, 
wanneer het systeem is ingebed in een warmtebad met temperatuur T. 


b. Beantwoord dezelfde vragen als het systeem door adiabatische wanden 
omgeven is. 


c‚ Wat voor verschil verwacht U als het gas niet ideaalis? 


12. We beschouwen een ideaal gas van N deeltjes met massa m in een -zwaarte- 
krachtveld g. De totale energie is de som van de kinetische energie K en 
de potentiele energie V. | 


waar z, de z-component van de positie van deeltje i is. Het gas heeft 
een temperatuur T, en bevindt zich in een cylinder met as evenwijdig 
aan de z-as, met grondoppervlak O op z = 0. 

Het bovenvlak van de cylinder bevindt zich op de hoogte h. 


a. Geef de toestandssom, de Helmholtz vrije energie, de warmtecapaciteit - 
bij een constant volume, de entropie en de chemische potentiaal in de 
limiet hee 


b, Wat is de gemiddelde hoogte van de deeltjes als ho 


c, Laat zien dat de dichtheid n als functie van z voldoet aan 


—Bmgz 


n(z) = n(o) e : 
Wegens de ideale gaswet geldt hetzelfde voor de druk 
—Bmgz | 
p(z) = p(o) e . 


(dit is de barometrische hoogteformule). 


d, Laat zien dat het drukverschil tussen twee hoogten z en z + dz juist over- 
eenkomt met het gewicht van de deeltjes tussen z en Z + dz. 

















13. a. Ga uit van de identiteit 
2 2 


ò F 9 F 
òVaT a TAV 


waar F de Helmholtz vrije energie is. 


Bewijs nu 


S d 
GP = GD 
T Vv 
b. Bewijs nu 
EE) En 


14. De Gibbs vrije energie is gedefinieerd als 
G = EE + pV - TS 
a. Bewijs nu uitgaande van de Emidanturete thermodynamische differentiaal 
dG = — SAT + Vdp + udN … 


Bij een homogeen systeem is de Gibbs vrije energie eenvoudig evenredig 
met het deeltjes aantal N, bij constante temperatuur en druk 

G(T,p, AN) = AG(T,‚p‚N) 
(Dit geldt bijv. niet voor een vat met gas en vloeistof samen.) 


b. Bewijs hieruit voor een homogeen systeem 


AG 
G(T‚p,N) = NGE 
T‚p 


c. Bewijs voor een homogeen systeem 
G(T,p,N) = iN .… 
1). a. Geef de groot-kanonieke toestandssom voor een klassiek ideaal gas. 


b. Bewijs dat voor een homogeen systeem geldt 
Y =spV ä 
waar Y de grand potentiaal is. 


c. Geef een vergelijking voor de druk van een klassiek ideaal gas, dat 
energie en deeltjes kan uitwisselen met een warmtebad. 


Î 

, 
là 
K 
8 
Ë 
î 

n 
Vis 
| 
EE: 
p 
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U WORDT VERZOCHT ELKE OPGAVE OP EEN APART VEL TE MAKEN 


Onderdeel Lorentz-invariantie 


ke 


Een elektron heeft een snelheid v = Cc. 


1 
Druk zijn totale energie E en ook de grootheid pc uit in de rustenergie E 


Lo 


| Ge mn ne 
Dit elektron maakt een | nt AS | foton 
kopbotsing met een foton (d.w.z. hun impulsen hebben tegengestelde 


richtingen). 


Bereken de energie die het foton minstens moet hebben (de drempelenergie) 
om de reaktie | Zan 


— + 
Yyte >e te + e@ 


mogelijk te maken. 
Geg.: E = 511 keV, 


In een gebied bevindt ‘zich een elektrisch geleidend medium waarin het — 
elektromagnetische veld en de ladings- en stroomdichtheden uniformen in 
de tija konstant zijn. 
X 
In het inertiale stelsel S (waarin het medium als geheel in rust is) is 
de stroomdichtheid In de ladingsdichtheid o, terwijl de geleiding voldoet 
aan de wet van Ohm | 
| er < > 
J = 0E 
je) Oo 


(o = de soortelijke geleiding). 
a) Voor een waarnemer in het inertiale stelsel S heeft het medium de snelheid u 


Wat zijn de stroomdichtheid J en de ladingsdichtheido in S? 


In stelsel S zal de wet van Ohm zeker niet luiden $ = 0 É. Geef hiervoor 


een kort argument. 


b) Als u = de h-snelheid van het medium en j = de h-stroomdichtheid in het 
medium beschouwen we s, gedefinieerd als | | 


S=} = u 
J en 


Les Ze 














Wat zijn de komponenten van s in S? 
LÄ LÀ | " LÀ, LÁ | s in SS? 
Bepaal 8 * u. 


Is s een tijdachtige of plaatsachtige h-vektor? 


c) Laat zien dat de vergelijking ie 4 | 
s=oFufe ad | | Dort: 
de relativistisch co-variante vorm van de wet van Ohm is | 
(F = benoe het É, B-veld). se 
Leid daartoe eerst af dat (1) in 5 inderdaad de wet van Ohm bevat, 


licht vervolgens kort toe waarom dit aart ook in andere enertiale stelsels 


ZO ÌS. 


Schrijf de wet van Ohm nu expliciet uit in S. 


Gegeven: transformatie van É als S' met snelheid v t.o.v. S beweegt: 


E' =ì E AE ne 
IM df | 


Veldtensor in korte notatie: 


(eÌ) - ik 
> 
1E O 
ji 


Onderdeel Statistische Fysika 


We beschouwen een systeem van N deeltjes, dat voldoet aan de Boltzmann 
verdeling. We nemen aan dat de toestandssom Z(N,V,T) voldoet aan de ge- 
lijkheid | 

8 ai 


Z(N,V,T) zr (Z,(V,T)) 


“waar 4 de toestandssom voor één deeltje is. RO geldt voor een ideaal gas). 


a) Geef een uitdrukking voor de Helmholtz vrije energie F en de chemische 
potentiaal u, in termen van Zee 
Gebruik de formule van Stirling 
log Ni = N log N = N, 


De grand potentiaal U is gedefinieerd door de gelijkheid. 


p= FE -— uN. 











ER ge OT DE EE EE EE EE mr ek NEER EN Winnend B OTE KE Ee 


'b) Laat zien dat 


y= — KIN 


Nu beschouwen we het systeem, terwijl het deeltjes kan uitwisselen met 


een reservoir, zodat het voldoet aan de Gibbsverdeling. 


c) Laat zien dat w voldoet aan dezelfde vergelijking, waar nu N het gemiddelde 


aantal deeltjes in het systeem is. 


Op het werkkollege is afgeleid 


Le = (®  , 
IV TN IT/V,N 


waar S de entropie, V het volume, p de druk en T de temperatuur is van een 
systeem met N deeltjes. U mag deze gelijkheid gebruiken zonder bewijs. 


Verder geldt de fundamentele thermodynamische differentiaal 
dE = TAS - pdV + uaN. 


a) Laat zien dat 


BE) ap 
Cv nor Pt+T Gm)yv 


Gebruik deze gelijkheden bij het volgende voorbeeld. Een van der Waalsgas 


van N deeltjes voldoet aan de vergelijking 


(p + a)(V — Nb) = NKT . 


1, 


Beschouw nu een van der Waalsgas dat met behulp van een wand opgesloten 
zit in de helft van een cylinder met volume V. Op zeker begintijdstip 
wordt de wand plotseling weggehaald en daarna zal het gas expanderen tot 
het totale volume V gevuld is. De cylinder is in warmtekontakt met een 


warmtebad van temperatuur T. 


b) Hoeveel is de aan het systeem toegevoerde warmte Q bij dit proces? 


zz EE 


ec) Hoeveel bedraagt de entropie-toename? 


